Determination of the single-ion anisotropy energy in a S" = 5/2 kagome 
antiferromagnet using x-ray absorption spectroscopy. 
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We report x-ray absorption and x-ray linear dichroism measurements at the Fe L2,3 edges of 
the geometrically frustrated systems of potassium and hydronium iron jarosite. Comparison with 
simulated spectra, involving ligand-field multiplet calculations modelling the 3d-2p hybridization 
between the iron ion and the oxygen ligands, has yielded accurate estimates for the ligand metal- 
ion hybridization and the resulting single-ion crystal field anisotropy energy. Using this method 
we provide an experimentally verified scenario for the appearance of a single-ion anisotropy in this 
nominally high-spin 3d^ orbital singlet ''S system, which accounts for features of the spin-wave 
dispersion in the long-range ordered ground state of potassium iron jarosite. 

PACS numbers: 75.25.-l-z, 71.70.Ch, 71.70.Ej, 75.10.Dg 



I. INTRODUCTION 

The jarosite group of minerals has received much at- 
tention for being close to ideal realizations of the kagome 
antiferromagnet. The kagome topology, which is named 
after a Japanese basket weaving pattern, frustrates the 
antiferromagnctic Nccl ordering^. Due to the resulting 
"underconstraint" neither in the quantum case oi S = 
1/2 nor in the classical limit is the kagome antiferromag- 
net expected to show a symmetry-breaking transition, 
even at T = K^i^i^. Physical realizations of this system 
are highly valued because they allow for the investiga- 
tion of the Mott insulating phase in the absence of Neel- 
like magnetic order. The jarosite group of general stoi- 
chiometry AM3(S04)2(OH)6, where A can be — amongst 
others — K, Na, Rb, Ag, NH3, or H3O, provides physical 
models for the kagome antiferromagnet with spin S = 
3/2 and 5/2, where the M site is occupied by trivalcnt 
Cr^i^i and Fe^i^iaiiaiii.ii^ii^ii^, respectively. M = V is 
a, S = \ system with ferromagnetic near-neighbor inter- 
actionsi^ii^. Fig. [1] illustrates the jarosite structure of 
well separated kagome layers consisting of MOe octahe- 
dra. Potassium iron jarosite is representative for most of 
the iron analogues. In this system a transition to a non- 
coUincar long-range ordered state is observed at 64 K, de- 
spite that the magnitude of the Curie- Weiss temperature, 
0CW, is as high as —800 K— li^. The ground state has 
shown to be a so-called g = state of positive chirality, in 
which the spins lie within the kagome plane and all point 
either in or out of each shared kagome triangle^ d^'^'^d^i^^ . 
The large figure 0cw/2n = 12.5 is an indication of rela- 
tively strong geometric frustration^^, but clearly the ob- 
served long-range ordered ground state points to addi- 
tional terms in the Hamiltonian of this kagome antifer- 



romagnet, beyond near-neighbor exchange. It has been 
shown that the q = Q ground state can arise due to mag- 
netic anisotropics such as easy-plane single-ion crystal 
field (CF) anisotropy DJl-E^y {Sl-Sl)^^w\ieT:e 
is the zero-field splitting parameter, or a Dzyaloshinsky- 
Moriya interaction (DMl)^^ D.jSi x Sj. The DMI term 
has shown to be symmetry allowed on the kagome lattice 
and has been argued to prevail over the CF anisotropy 
in the spherical (L = 0) 3d^ Fe^+ ior^. Fits to the 
spin-wave dispersion curves measured with neutron spec- 
troscopy on single crystal^^ and powder samples^ of 
potassium jarosite have found to be slightly better in 
a model with DMI compared to the model with CF 
anisotropy. The evidence so far has, however, been in- 
conclusive. 

Given the similarity in crystal structure of potassium 
(iron) jarosite and hydronium (iron) jarosite the above 
effects could expected to be equally important in the lat- 
ter. However, hydronium jarosite has a spin-glass ground 
state with a freezing temperature of Tg = 17(2) K2^. 
The Weiss temperature for hydronium jarosite is less 
accurately known (—700 to —1400 K) because the in- 
verse susceptibility does not enter a linear regime suffi- 
ciently far below the decomposition temperature of the 
compoundiS. It has been suggested that the spin-glass 
state is the result of proton transfer from the hydronium 
groups to the Fe-(OeqH)-Fe super-exchange mediating 
hydroxy groups^. One indication that this could happen 
might be that the OH groups at the crystallographically 
analogous location in the related kagome compound zinc 
paratacamite can form a muonium state by chemically 
binding a positive muon^^. The aim of the work de- 
scribed here is to investigate the origins of the magnetic 
anisotropy (of either single-ion or DMI character) in the 
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iron jarosites, and to understand the essential differences 
between tlie salts such as potassium jarosite with a long 
range ordered ground state and hydronium jarosite. 




FIG. 1: (color online) The jarosite structure viewed perpen- 
dicular to the c axis (left picture) and along the c axis (right 
picture), revealing the kagome network of magnetic ions. The 
MOe octahedra are drawn in transparent brown with the 
M^^ ions at the center in dark brown. The SO4 tetrahe- 
dra are transparent yellow. The large spheres at the origin of 
the unit cells in blue are the A site ions, and the hydrogen of 
the OH groups are drawn as small white spheres. 

Using Mossbauer spectroscopy22i2S the spin state of the 
Fe'^"'" in the iron jarosites was found to be 5' = 5/2, thus 
nominally high-spin 3d^. The free-ion state is ^S, 
i.e. L = 0. Hence, there is no orbital angular momen- 
tum due to the spin-orbit coupling in the presence of a 
crystal field as long as the system is high-spin 3dP. This 
holds at least up to third order perturbation theory of 
the trigonal crystal-field distortions, spin-spin and spin- 
orbit couplings^. Even in a crystal field the 3dP 
S = 5/2 spin should therefore be practically isotropic 
apart from the possible effects of dipolc-dipole interac- 
tions. Both DMI and single-ion anisotropy arise only 
when the spin-orbit operator L ■ S can induce a finite 
orbital angular momentu m^^i'^" . That (L) > in the 
iron jarosites is clear from Curie- Weiss fits to the high 
temperature magnetic susceptibility which indicate 
between 6.2 and 6.7 in potassium iron jarosite^ and 
6.6(2) /iB in the hydronium analogu o^^i'^^ (compared to 
a spin only value of the effective moment of 5.92 fiB)- 
Clearly, this opens up the possibility for both DMI and 
single ion anisotropics. However, the origin of the or- 
bital angular momentum itself has so far received little 
attention. 

A very strong crystal field, such as in iron phthalocya- 
nina^ can give rise to a low-spin ground state in which 
the spin-orbit coupling will in general reinstate part of 
the orbital angular momentum of the free-ion configura- 



tion which is, in zero order, quenched by the crystal field. 
It could be that in weaker crystal fields, such as those 
present in iron oxides, some orbital angular momentum 
appears due to a mixing in of the low-spin state into the 
ground state. The crystal field alone does, however, not 
provide matrix elements between high-spin and low-spin 
states. The most likely explanation is therefore that an 
Fe^"*" configuration mixes in due to charge transfer— from 
the surrounding ligands. For hematite, e.g., the effective 
valence of the Fe ion has been estimated at 2.43 due to 
charge transfer from the oxygen ligands^. These effects 
have earlier been shown to explain zero field splitting in 
Mn2+ and Fe^+ substituted crystals^i^^ and in the Mn^+ 
compounds MnF^S, MnPSa and MnPSeg^, despite the 
lower covalcncc of the fluoride and phosphorus-containing 
ligands compared to oxygen. 

Transition metal 2p — > 3d x-ray spectroscopy studies 
have played a key role in the understanding of the ef- 
fect of charge-transfer and crystal-field in determining 
the electronic and magnetic properties in transition met- 
als and their compounds^'^^i^'^. We have measured the 
x-ray absorption and x-ray linear dichroism (XLD) at the 
Fe L2,3 (2p — > 3c?) edges on single-crystalline and pow- 
der samples of potassium and hydronium iron jarosite. 
By comparison with atomic multiplet calculations, tak- 
ing into account ligand-ficld effects and multiple ionic 
configurations, we obtain accurate values for the charge 
transfer and single-ion anisotropy in the 5 = 5/2 kagome 
antiferromagnet. The results are compared to magnetic 
susceptibility and neutron spectroscopy measurements on 
in particular potassium iron jarosite. 



II. EXPERIMENTAL 

Single crystals of potassium iron jarosite, 
KFe3(OH)6(S04)2, were grown using a hydrother- 
mal reduction-oxidation method as described in Ref. 
4.88 g (28.0 mmol) K2SO4 and 2.2 niL (40 mmol) H2SO4 
were dissolved in 50 mL distilled water, and transferred 
into a 125 mL PTFE liner of a stainless steel bomb. 
0.56 g (10 mmol) iron wire with a diameter of 2 mm 
was added to the solution. The bomb was placed in an 
oven at 202°C for 4 days, then cooled down to room 
temperature at a rate of 0.3°C/min. The precipitate was 
washed, filtered and dried, yielding 0.37 g, which is 22% 
based on Fe. In order to obtain sufficiently large crystals 
it was important that the surface area of the solution 
exposed to air inside the vessel was minimized. This was 
achieved using a glass container which fitted inside the 
PTFE liner, with only a small hole at the top. Using this 
method single crystals of up to 0.7 mm in diameter were 
obtained, which was sufficiently large for the application 
in polarized x-ray spectroscopy. A number of crystals 
were characterized using single-crystal x-ray diffraction. 
These crystals were identified as potassium iron jarosite 
and face-indexed. 

The original solvothcrmal synthesis method^i^ 
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for the jarositcs was used for the preparation of 
H30Fe3(OH)g (804)2- In this preparation 6.6 g 
(22 mmol) of Fe2(S04)3-5H20 was dissolved in 50 mL 
water. This solution was transferred to a 125 mL PTFE 
liner of a stainless steel bomb. The solution was heated 
in the bomb to 140°C for 12 hours. The hydronium 
iron jarosite which precipitated during the reaction was 
washed, filtered and dried, yielding ^^0.27 g of product. 

Polarized soft x-ray absorption measurements were 
carried out using the liquid helium cryostat, housed in the 
high-field superconducting magnet on ID08 at the Euro- 
pean Synchrotron Radiation Facility (ESRF) in Greno- 
ble, France. Total-electron-yield spectra were measured 
by recording the drain current from the sample as a func- 
tion of photon energy in the region of the Fe L2,3 absorp- 
tion edges. Potassium iron jarosite is strongly insulating, 
so the resonance enhanced total-electron-yield signal at 
the iron edge is much smaller than the signal from the 
surrounding metal of the sample holder. To avoid this 
source of background intensity a potassium iron jarosite 
single crystal of 0.7 mm diameter was attached using 
silver glue to the end of a graphite tip, which in turn 
was glued to the end of a standard ESRF sample holder. 
The sample was aligned in the x-ray beam with the crys- 
tallographic c axis pointing in the horizontal plane and 
perpendicular to the incident beam. This allowed us to 
measure absorption spectra with the incident x-ray po- 
larization both perpendicular and parallel to the c axis. 
Spectra were recorded with alternating horizontal and 
vertical polarization, and for each polarization averages 
were taken over at least eight scans. In this way the 
polarization dependent spectra were taken at a number 
of temperatures between 20 and 290 K, in zero external 
field. 

Isotropic absorption spectra at the Fe L2.3 edges were 
obtained on beamline 5U.1 of the Synchrotron Radiation 
Source (SRS) at Daresbury Laboratory, UK, using potas- 
sium and hydronium iron jarosite powder dispersed over 
double-sided UHV-compatible carbon tape. 



transition giving rise to a I'Pxji core hole (Li edge) is 
at ~12 eV higher photon energy due to the 2p spin-orbit 
couphng. The inset of Fig. [5] shows the /c and /ab spectra 
in close up at the L3 edge, which is the most significant 
region of the spectrum. 
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FIG. 2: (color online) x-ray absorption spectra of potassium 
iron jarosite at 290 K measured with x-ray polarization par- 
allel and perpendicular to the c axis, 7c (black circles) and 7ab 
(red crosses), respectively. The lower trace shows the corre- 
sponding values of the isotropic absorption, Jiso = f^ab + |7c 
(black line) and the x-ray linear dichroism, Ixm = 7ab — 7c 
(thin black line) . Also shown 7iso from a powder sample mea- 
sured on 5U.1 at SRS Daresbury with a lower energy resolu- 
tion (broken, red line). The inset shows a close up of 7c and 
7ab at the L3 



III. RESULTS 

Figure[2]shows the spectra obtained from single-crystal 
and powder samples of potassium iron jarosite at 290 K. 
The upper curves in Fig. [2] display the x-ray absorption 
spectra from a single crystal measured on beamline ID08 
at the ESRF, with the incident x-ray polarization parallel 
to the c axis, Ic (black circles), and perpendicular to 
the c axis, /ab (red crosses). The curves in the middle 
of the figure give the corresponding isotropic absorption 
spectrum /jgo = f ^ab + §-^c (black line) and the isotropic 
spectrum directly measured from a powder sample on 
station 5U.1 at the SRS (broken red line). The curve at 
the bottom gives the x-ray linear dichroism /xm defined 
as /ab — /c (thin black line). The 2p 3d transition 
consists of two edges, the first one at ~709 eV (L3 edge) 
corresponds to the creation of a 2p3/2 core hole. The 



As seen in Fig. [31 the structure at the L3 edge be- 
comes more prominent in the isotropic spectrum as the 
temperature is lowered. At each stabilized temperature, 
the reproducibility of up to eight spectra demonstrated 
the absence of any temporal variation in the spectral line 
shape or intensity, therefore charging effects and drift in 
sample alignment can be eliminated as the cause of the 
observed temperature dependence. These isotropic spec- 
tra obtained by weighted averaging of /ab and are in 
good agreement with the spectra from powder samples 
measured for the SRS at all temperatures. That this 
change in line shape occurs in the powder averaged spec- 
trum might be an indication of a slight change in the 
electronic structure as the temperature is lowered. Only 
very small changes were observed in the XLD spectral 
line shape on cooling the sample, as shown in Fig. 31 
This change in the spectral line shape is accompanied by 
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a slight reduction in the integrated XLD spectrurc4i 

4id(w)dw/ / /iso(w)dt^ (1) 



corresponding to a small change of the total quadrupole 
moment of 0.05(2). This might be interpreted as that the 



Fe 3d shell gains between 3 and 7% net 
at base temperature. 
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FIG. 4: (color online) Temperature dependence of the x-ray 
linear dichroism measured at the Fe L2,3 edges from a syn- 
thesized potassium iron jarosite single crystal (black lines), 
compared to calculated linear dichroism corresponding to the 
isotropic simulated spectra of Fig. |3l scaled in intensity to fit 
the experimental spectra by factors ranging from 1/2 at 290 K 
down to 1/3 for the 60 and 40 K data. Possible reasons for 
this mismatch are discussed in the text. 



FIG. 3: (color online) Temperature dependence of the 
isotropic x-ray absorption measured at the Fe L2.3 edges from 
a synthesized potassium iron jarosite single crystal (black 
lines), compared to calculated spectra (as discussed in Sees. 
IV and V) with simulated temperatures as listed in the legend. 



The room temperature spectra of hydronium iron 
jarosite were identical to the spectra obtained on potas- 
sium iron jarosite within the experimental resolution as 
measured at the SRS (the red, broken line in the main 
panel of Fig. [2] for comparison). 



IV. LIGAND-FIELD MULTIPLEX 
CALCULATIONS 

Due to the strong electrostatic interactions between 
the 2p core hole and the 3d levels in the final state, the 
2p — > 3d absorption spectrum is not simply proportional 
to the density of unoccupied 3d levels^ as a function of 
energy. For ferromagncts and magnetically soft materi- 
als, where all the magnetic moments can be aligned by 



applying an external field, the x-ray magnetic circular 
dichroism (XMCD) and x-ray magnetic linear dichroism 
(XMLD) sum rules can be applied to obtain the expecta- 
tion values of the orbital and spin magnetic moments^ 
and the magneto-crystalline anisotropy^, respectively. 
In the case of the jarosites it is not possible to align the 
magnetic moments parallel with an external field due to 
the large energy scale of the antiferromagnetic interac- 
tions (6cw ~ —800 K). The only way to obtain accurate 
information about the relevant expectation values is by 
explicitly calculating the transition probabilities 



(2) 



between the ground state ipg and excited state i/'e , where 
fq is the electric-dipole operator for x-ray absorption 
with polarization q, and F is the lifetime broadening, 
which is treated as a fitting parameter. The experi- 
mental x-ray spectra were used to obtain good approx- 
imations for tpg and ipe with the aid of the Hilbert++ 
code^i^. The Hilbert space for \ipg) is spanned by deter- 
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minants of the lowest atomic LSJ configurations of |3d") 
and |3(i""'"^L), where L denotes a hgand hole. Likewise, 
the Hilbert space for \'ijje) is spanned by determinants of 
the atomic configurations of |2p^3(i"+^) and |2p^3d"+^i). 
The Slater integrals for the Zd-id and 2p-3d interac- 
tions and the 2p and 3d spin-orbit parameters for the 
ground state and final state configurations calculated us- 
ing Cowan's atomic multiplct program^ are tabulated in 
Ref. [U. Approximate solutions for the eigenfunctions of 
the iron ion in the compound are then obtained by diag- 
onalization of the Hamiltonian including charge-transfer 
and crystal-field terms arising from the surrounding oxy- 
gen ligands^ 

^ "^atom 

+ Yl^t^^bd+^'-r^Pi + ^^-bidtzPi + <4zPy) + Cc] 
b 

b 

+ <^pJ2^PtPi +PtPy +PtPi)- (3) 
b 

In this Hamiltonian, the first term represents the atomic 
Hamiltonian. The second term gives the hybridization 
between metal 3d and oxygen 2p orbitals, where t^-.b 
and <7r,& are the Slater-Koster hopping parameters. For 
each bond, ta-^T^ are rescaled with the bond length using 
ta,Tv,b = ta,TT{Rh/ Rrci)°' whcrc Ri, / R-^cf is the normalized 
length of Fe-0 bond b and the rescaling exponent a is 
treated as a free parameter. The third term gives the 
contribution of the electrostatic crystal field which was 
not used in our case. The last term gives the energy of 
the 2p electrons in the oxygen valence band, where €p 
is a simulation parameter related to the charge-transfer 
gap Apd via Ap^ = U^d^ — tp. The total dd Coulomb 
energy for the addition of a sixth electron in the Fe id 
shell C/sds = Y^n=i '^^^ found to be 44.3 eV by con- 
sidering the case tc^-n = 0, where no charge transfer oc- 
curs unless the Fe^"*" d^L configuration has a lower energy 
than the Fe'^^" d^ configuration. The sum is made over 
all bonds 6, where the local coordinate frame x,y,z is 
oriented with the z axis along the bond direction. Nu- 
merous comparisons between experiment and theor y^^i^^ 
have shown that this approach works well, provided that 
the Slater integrals for the Coulomb interactions and 
spin-orbit couplings are reduced to 70-80% of the atomic 
Hartree-Fock values^. Here, all pd and dd Slater inte- 
grals were reduced to 80% of the atomic values calculated 
using Cowan's program. The lifetime broadening F was 
set to 0.24 eV for the L3 edge and 0.3 eV for the L2 edge. 
These are values as commonly used in the calculation of 
x-ray absorption spectra, with L2 always slightly larger 
than L3. This spectral broadening is due to the coupling 
with the (infinite dimensional) electromagnetic field. 

The crystal field at the iron atom in potassium iron 
jarosite is approximately of D^h symmetry, with an elon- 
gation of the bond length i?ap along the fourfold symme- 
try axis (c') of ^^4% with respect to the bond length i?cq 
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FIG. 5: (color online) Experimentally obtained isotropic spec- 
trum and x-ray linear dichroism (black lines) at 290 K and 
simulated XAS and XLD (red lines) Boltzmann averaged 
at 110 K. The simulated XLD is roughly twice as large as 
the experimentally obtained XLD which is ascribed to mis- 
alignments and sample imperfections. 



in the equatorial planers. The local symmetry axes of the 
FeOe octahedra are canted with respect to the crystallo- 
graphic c axis by an angle of 19.2°, resulting in buckled 
Fe04 planes forming the kagome layers. 

Our simulation shows that the spectra are mainly de- 
termined by the strongly anisotropic cr-type hybridiza- 
tion between the iron 3d and oxygen 2p orbitals. The 
electrostatic part of the crystal field could be left out. 
An excellent agreement between experimental and sim- 
ulated spectra for Iiso was obtained using only four 
simulation parameters, namely the charge-transfer gap 
Apd; ^(T, and the rescaling constant a as defined in 

CP, '= t°<^,(i?ap/i?eq)" (with i?ap/i?oq = L035 in the 

present case). When the simulated isotropic spectra is in 
good agreement with the experimental spectrum, at the 
same time optimal agreement is obtained for the linear 
dichroism spectrum. At high temperature there is good 
agreement between the experimental and simulated lin- 
ear dichroism spectrum, apart from an intensity scaling 
factor of 0.5 as shown in Figs. [4] and [5] With the fitting 
parameters for the simulation there is no way to reduce 
the overall intensity of the linear dichroism without rad- 
ically modifying the isotropic spectrum. For up to 6% 
the intensity difference might be attributable to the local 
canting of the FeOe octahedra with respect to the ah- 
plane. In addition to this there is an inevitable loss of 
intensity in the XLD due to misalignments and sample 
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imperfection. Reduced energy resolution can also lead 
to a reduction of the measured XLD. We conclude that 
the difference in overall intensity of the experimental and 
calculated /xm is not due to fundamental shortcomings of 
the simulation. For reasons we do not quite understand, 
the agreement between theory and experiment as evident 
in the temperature dependence of the isotropic spectra 
is not reflected in the XLD spectra at lower tempera- 
tures. We expect that these spectra, which are averages 
over increasingly fewer states, are also increasingly sensi- 
tive to shortcomings in our simulations. Since the XLD 
is an orientational energy dependence it will be partic- 
ularly sensitive for any shortcomings in our theoretical 
model. It should be emphasised that, despite the differ- 
ences between the experimental and simulated spectra 
at low temperatures, the overall agreement between ex- 
periment and simulation as obtained here is remarkably 
good. The parameters as found in our simulation are 
given in Table [T] (left-hand side). The large error bar on 
the value for the charge-transfer gap of 3.5(2.0) eV re- 
flects the adjustment of the total dd Coulomb energy on 
the iron site ^7^^' — {/^^^{{n) — 5) to any increase in Cp 
by transfer of additional electron weight into the Fe ion. 
This is reflected in the total error in (n) of '^0.1 electron. 

The values found here for Apd and are similar to 
those found previously for other Fe'^+ compounds^i^. 
FC2O3, which is in many respects probably the com- 
pound most closely related to iron jarosite, was previ- 
ously found to have a band gap of predominantly charge- 
transfer character of 2 to 2.5 eV— . With an unusually 
small value for of 0.3(2) eV the lowest energy peak, 
at 707.8 eV is reproduced in the simulated spectrum. 
The main peak at 709 eV splits into two peaks with the 
introduction of a trigonal distortion of the FeOg octahe- 
dron. A common value for the rcscaling exponent a is 
—3 and the small value found here (—7.5) is an indication 
that the charge anisotropy is larger than the crystal field 
anisotropy as defined by the oxygen positions around the 
iron. The lowest energy low-spin configuration, corre- 
sponding to the *Fii/2 free-ion configuration, lies 1.2 eV 
higher in energy and is not found to mix into the ground 
state. 

The singlc-ion anisotropy Dz is obtained as the dif- 
ference between the energy of the Sz = ±1/2 and ±5/2 
levels and was found to be of the easy-plane type. This 
easy-plane anisotropy was also confirmed by rotating the 
spin quantisation axis to various directions with respect 
to the crystal field. In this way we have also attempted 
to obtain estimates for E^y The D^h symmetry of the 
ligand field was lowered to C2 , to model the actual posi- 
tions of the oxygen anions around the cation as obtained 
from crystallographical datapi^ii^. However, this had no 
effect on the agreement of the calculated spectra with ex- 
periment, nor on the values as listed in Tableland hence 
we consider an estimate for E^y from our data presently 
not attainable. It should in this respect be noted that 
the value for was found to change with the splitting 
of the main peak at 709 meV, introduced by the trigo- 



TABLE I: Results from simulations of the Fe 1/2,3 spectra 
of potassium iron jarosite, with ligand-field multiplet calcu- 
lations using the Hilbert++ code^. The left-hand side lists 
all the parameters varied in the simulation. The right-hand 
side gives the resulting expectation values for the calculated 
ground-state properties of the Fe ion. The errors indicate the 
range over which the parameters could be varied whilst main- 
taining a good agreement with the data. The right-hand side 
of the lowest panel gives the effective moment in the case of 
full orbital polarization, while the calculation indicates only 
a very weak orbital polarization of the added electron weight 
due to charge transfer. We use the convention that upper case 
J, L, S represent operators and lower case j, I, s represent the 
quantum numbers given by (L^) = l{l + 1) etc. 
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nal distortion of the ligand field. Its value is therefore 
tightly constrained by the experimentally obtained spec- 
tra and the 3d spin-orbit coupling strength. As is detailed 
in Table [J our calculation indicates that for the isolated 
FcOg cluster most of the orbital angular momentum of 
the Fe 3d^ configuration is quenched by the ligand field, 
with (L) = 0.008(3) while the maximum orbital angu- 
lar momentum which could in this case be realised as 
7(1^=1.5(2). 

Due to the lower resolution at which the hydronium 
jarosite powder spectrum has been measured we cannot 
provide detailed information about this system, apart 
from a confirmation that the general crystal field is sim- 
ilar, with comparable charge-transfer and (cr-type) hy- 
bridization. This is a good indication that the magnetic 
exchange in hydronium jarosite is likely to be comparable 
with that of the other jarosites, i.e. with a Weiss temper- 
ature around —800 K. There is also a slight evidence that 
the ligand-field anisotropy and hence the zero-field 
splitting, is slightly smaller in hydronium jarosite. 



V. DISCUSSION 

The expectation values of the ground-state obtained 
in our calculation are given in Table [J (right-hand side). 
The corresponding free-ion configuration is d'^ with 
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an additional electron weight of '^0.4 electron that is 
mainly in the d{x^ — y^) orbital. The spin-orbit cou- 
pling is much smaller than the t^r hopping bandwidth 
and hence the fraction of the electron transferred to the 
ligands to the Fe cation is only very weakly orbitally po- 
larized, with {L) = 0.008(3). This value is much smaller 
than the orbital angular momentum as inferred from 
high-tcmpcrature magnetic susceptibility data, which in- 
dicate that fieS ~ 6.5 iig^. As shown in the lower panel 
of Table [I] (right-hand side) the experimental value is 
remarkably close to the case where the added electron 
weight due to charge transfer is fully orbitally polarized. 
It is not expected that in this 3d transition metal com- 
pound I and j{~ s + /) are good quantum numbers, but 
there is sufficient evidence that (L) is significantly larger 
than zero^. It might be that comparison of experimen- 
tal results with the Mott-Hubbard-likc calculations as de- 
scribed here, but with larger clusters containing two or 
three iron centres, can provide an explanation for the 
discrepancy between the calculated and experimentally 
obtained orbital angular momentum. 

The energy of = 0.5(1) meV for the easy-plane 
single-ion anisotropy is in good agreement with the values 
found by Matan et al^ (0.428(5) meV) and Coomer et 
al^ (0.47(2) meV) in their spin- wave analysis of inelastic 
neutron data from single crystals and powders, respec- 
tively. This implies that an easy-plane anisotropy, along 
with further neighbor interactions, can explain the mag- 
netic ground state in potassium iron jarosite. Further ev- 
idence might be found in the temperature dependence of 
the experimental and calculated isotropic spectra. The 
spectra in Fig. [5] were taken at room temperature and 
were compared with the calculated spectra for the Boltz- 
mann averaged configurations of the FeOg cluster. As is 
visible in Fig. [3] the relative intensity of the two maxima 
in the split peak at 709 eV changes in a similar way as 
was observed in the experimentally obtained spectra. At 
the same time an (unexplained) increase in intensity of 
the peak at 707.8 eV peak relative to the intensity of the 
main peak is observed in the experimental spectra. We 
have not been able to account for the latter (which was 
also observed in powder samples) either by changing the 
temperature of the calculated spectra or with simulations 
with new values for Spd, io-.Tr and a. The change in the 
main peak at 709 eV where we obtain good agreement, 
is a direct result of the spin-orbit coupling which causes 
a change in the shape of the electronic 3d orbital as the 
spins align in the easy plane. The change in shape of the 
main peak in the experimental spectra was found to occur 
between 200 K and 40 K. This is in agreement with pre- 
vious experiments which have shown that in potassium 
jarosite the alignment of spin with the kagomc planes sets 
in above 120 K, i.e. well above the transition temperature 
of 65 KH. The temperature dependence in the experi- 
mentally obtained spectra was simulated with the tem- 
perature (given in units of eV) used for the Boltzmann 
averaging over excited states as the only extra parame- 
ter, as shown in Fig. [31 All other simulation parameters 



were fixed to the values given in Table HI In this way 
we found the expectation values for (|S'z|) at each cor- 
responding experimental temperature, arising from the 
thermal population of the ms levels in the 5 = 5/2 mul- 
tiplet. The temperature dependence of {\Sz\) (Fig. [6]) 
shows a gradual alignment of the spins into the kagome 
planes, at temperatures well above the transition tem- 
perature to a long-range ordered state at 64 K. This is 
in rough agreement with the degree of co-planarity of 
the spins as measured using neutron diffraction on large 
single crystalsi^. 

2| 1 
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FIG. 6: The expectation values {\Sz\} of the simulated spec- 
tra, against the temperature of the corresponding experimen- 
tal spectrum. For a single S — 5/2 spin the minimum for 
(ISzl) is 0.5 and the maximum is 1.5 (equal population of all 
m$. The gradual increase of (IS^I) with temperature is in 
agreement with earlier neutron measurements which show a 
gradual alignment of the spins in the kagome plane at tem- 
peratures well above the transition temperature of 64 K^. 



The simulated linear dichroism spectra at the low- 
est temperatures, corresponding to states with the spins 
aligned in the equatorial plane (the states Is,™^') = 
|2.5, ±0.5)) do not agree well the experimentally obtained 
spectra as shown in Fig. IH although some features in 
the temperature dependence of the XLD are reproduced 
in the simulated spectra. Without success we have at- 
tempted to improve on this situation by calculating the 
spectra for the Fe spin aligned along the x' axis in the 
equatorial plane (i.e. |s, m^') = |2.5, ±2.5)) which resem- 
bles more closely the actual ground state in potassium 
jarosite. The remaining disagreement could be in part 
due to experimental error due to sample imperfections 
and in part due to intrinsic limitations of the simulation. 

It might be naive to expect that the ligand-field can be 
inferred from crystallographic data alone. A first indica- 
tion for this is the large value of a obtained here, which 
amplifies the effect of the crystallographic trigonal dis- 
tortion of the FeOg octahcdra. Even when the symmetry 
of the ligand field is further lowered from D^h and accu- 
rately modelled on the locations of the oxygen anions 
around the central iron^ neither the calculated spec- 
tra nor the simulation parameters change significantly. 
Therefore, it is possible that the actual charge anisotropy 
follows more closely the xy plane rather than the equa- 
torial oxygens. The easy-plane magnetic anisotropy of 
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= 0.5(1) meV found here (see Table is then not 
directly related to the trigonal distortion of the crystal 
field but more to the 2D character of the electronic hop- 
ping and magnetic exchange, the virtual absence of tt 
symmetry hybridization and the Fe 3d spin-orbit cou- 
pling strength. In that case the easy plane could coin- 
cide more closely with the kagomc plane. The strong 
2D character can in part be attributed to the sulphate 
groups separating the kagome layers. The formal valence 
of the sulphur cation is 6-1- with the [Ne] electronic con- 
figuration, but this large electrostatic charge is likely to 
attract electron weight from the neighboring oxygen an- 
ions, including the oxygens apical in the FeOg octahedra. 
The hopping parameters ta.-R from these oxygens to the 
iron is then further reduced with respect to the in-plane 
hopping. This translates into the large value for a found 
here. 

The large difference between and obtained here 
might be explained by the triangular arrangement of the 
neighboring FeOe octahedra. This dramatically lowers 
the symmetry of the local crystal field around the oxy- 
gen ligands and reduces the tt symmetry overlap of the 
oxygen 2p orbitals with the iron Sd orbitals. Due to the 
low symmetry of the oxygen sites the symmetry of the 
crystal field can be expected to be even lower than indi- 
cated by the crystallographic positions of the anions. It is 
in this sense fortunate that the situation can be modelled 
so well by the anisotropic tt and a hopping whose ratio is 
equal for all six anions around the central iron. It could 
well be that the last differences between the calculated 
and experimental spectra, such as the low temperature 
x-ray linear dichroism, some of the difference in the over- 
all intensity of the XLD and the missing bump just above 
707 eV, could be brought into agreement by allowing a 
different balance between and 1^ for the apical oxygen 
anions. 

In neutron spin-wave studics^ii^i^ the crystal field 
anisotropy scenario has been compared in close detail 
with that of DMI. In the present context such a discussion 
would require calculations on clusters containing two to 
three iron centres. Such a calculation would certainly be 
of great interest, but has not been attempted here, since 
the present approach already accounts so well for the ex- 
perimental observations on potassium iron jarosite. The 
present result closes the last missing causal link, from 
electronic structure to collective magnetic ground state. 

VI. CONCLUSION 

The isotropic and x-ray linear dichroism spectra taken 
at the Fe L2,3 edges in potassium iron jarosite indicate 
that an orbital angular momentum and zero-field split- 
ting arise due to an anisotropic, mainly a type, charge- 



transfer from the oxygen ligands. This closely resembles 
the scenario described in Ref. [H to explain the zero-field 
splitting as observed in Mn^+ salts. 

The zero-field splitting parameter of the crystal field 
anisotropy of 0.5 meV as found here is in excellent 
agreement with the value found from fits of a crystal 
field anisotropy model to the spin-wave spectrum in the 
ground state of potassium iron jarosite2ii22iiH. However, 
the most important result obtained here is that we have 
verified in detail the mechanism by which this magnetic 
anisotropy arises, closing the remaining missing link be- 
tween the electronic structure in potassium iron jarosite 
and its magnetic ground state. We have even mea- 
sured a change in the electronic structure as the mag- 
netic ground state becomes co-planar, with the change 
of the isotropic line shape over the relevant temperature 
range. Of course, an additional in-plane anisotropy E^y 
and further- neighbor interactions are needed to fully 
explain the magnetic ground state. The present study 
does not provide more hints to the origin of the glassy 
state in hydronium jarosite. We found that the energy 
scale of the magnetic interaction in hydronium jarosite 
should be comparable to that in potassium jarosite. It 
might be that temperature dependent spectra taken at a 
higher resolution could provide new insights in the case 
of hydronium jarosite. A final interesting observation is 
that the maximum orbital angular momentum polariza- 
tion possible for the amount of charge-transfer observed 
here, but not the actually calculated orbital polarization, 
is in rough agreement with the experimentally observed 
effective moment. 
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